On the mechanism for orbital-ordering in KCuFa 
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The Mott insulating perovskite KCuFs is considered the archetype of an orbit ally-ordered system. 
By using the LDA+dynamical mean- field theory (DMFT) method, we investigate the mechanism for 
orbital-ordering (00) in this material. We show that the purely electronic Kugel-Khomskii super- 
exchange mechanism (KK) alone leads to a remarkably large transition temperature of Tkk ^ 350 K. 
However, orbital-order is experimentally believed to persist to at least 800 K. Thus Jahn- Teller 
distortions are essential for stabilizing orbital-order at such high temperatures. 

PACS numbers: 71.10.Fd, 71.10.-w,71.27.+a,71.10.Hf 



In a seminal work [l[ Kugel and Khomskii showed 
that in strongly correlated systems with orbital degrees 
of freedom many-body effects could give rise to orbital- 
order (00) via a purely electronic super-exchange mech- 
anism. Orbital-ordering phenomena are now believed to 
play a crucial role in determining the electronic and mag- 
netic properties of many transition-metal oxide Mott in- 
sulators. While it is clear that Coulomb repulsion is a 
key ingredient, it remains uncertain whether it just en- 
hances the effects of lattice distortions [2| or really drives 
orbital-order via superexchange [i^. 

We analyze these two scenarios for the archetype of 
an orbitally-ordered material, KCuFa [1]. In this ZdP 
perovskite the Cu (i-levels are split into completely filled 
three-fold degenerate t2^-levels and two-fold degenerate 
e^-levels, occupied by one hole. In the first scenario 
Jahn- Teller elongations of some Cu-F bonds split the par- 
tially occupied e^-levels further into two non-degenerate 
crystal-field orbitals. The Coulomb repulsion, [/, then 
suppresses quantum orbital fluctuations favoring the oc- 
cupation of the lower energy state, as it happens in some 
t2^-perovskites [3|, (j]. In this picture the ordering is 
caused by electron-phonon coupling; Coulomb repulsion 
just enhances the orbital polarization due to the crystal- 
field splitting [sl, Isf. In the second scenario the purely 
electronic super-exchange mechanism, arising from the 
e^- degeneracy, drives orbital-ordering, and Jahn- Teller 
distortions are merely a secondary effect. In this picture 
electron-phonon coupling is of minor importance [1]. 

The key role of Coulomb repulsion is evident from 
static mean-field LDA+/7 calculations, which show 0,0] 
that in KCuFa the distortions of the octahedra are sta- 
ble with a energy gain ~ 150 — 200 meV per for- 
mula unit, at least an order of magnitude larger than 
in LDA [e", ^ and GGA 0, 0; recent GGA+DMFT 
0] calculations yield very similar results, suggesting in 
addition that dynamical fiuctuations play a small role 
in determining the stable crystal structure of this sys- 
tem. However, these results might merely indicate that 
the electron-phonon coupling is underestimated in LDA 
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FIG. 1: (Color online) Left: Crystal structure and orbital- 
order in a- type 12] KCuFs. Cu is at the center of F octhaedra 
enclosed in a K cage. The conventional cell is tetragonal with 
axes a, b, c, where a — h^ c = 0.95aV2- The pseudo-cubic axes 
are defined as x = (a+b)/2, y = (— a+b)/2, and z = c/2. All 
Cu sites are equivalent. For sites 1 the long (short) bond / (s) 
is along y (x). Vice versa for sites 2. Orbital |2) (see table |l|, 
occupied by one hole, is shown for each site. Right: Jahn- 
Teller distortions at sites 1, measured by ^ = (/ — s)/(/ + s)/2 
and 7 = c/a\/2. R is the experimental structure, R5 and I5 
two ideal structures with reduced distortions, and /o is cubic. 



or GGA, probably due to self-interaction, rather than 
identifying Kugel-Khomskii super-exchange as the driv- 
ing mechanism for orbital-order. This is supported by ab- 
initio Hartree-Fock (HF) calculations which give results 
akin to LDA+/7 [9]. Moreover, in the superexchange sce- 
nario it remains to be explained why Too ^ 800 K [lo[, 
more than twenty times the 3D antiferromagnetic (AFM) 
critical temperature, Tn ~ 38 K [HI [13, a surprising fact 
if magnetic- and orbital-order were driven by the same 
super-exchange mechanism. 



2 



In this Letter we study the Kugel-Khomskii mechanism 
at finite temperature and identify the origin of orbital- 
ordering in KCuFs. We will show that super-exchange 
alone leads to orbital-order with Tkk ^ 350 K, less than 
half the experimental value. Thus Jahn- Teller distortions 
are essential for driving orbital-order above 350 K. 

KCuFa is a tetragonal perovskite made of Jahn- Teller 
distorted C uFp octahedra enclosed in an almost cubic 
K cage [m, [3- The Jahn- Teller distortion amounts to 
a 3.1% elongation/shortening of the CuF distances in 
the xy-plane. The tetragonal distortion reduces the CuF 
bond along z by 2.5%, leaving it of intermediate length. 
The long (/) and short (s) bond alternate between x and 
y along all three cubic axes (a- type pattern) [l2| . At each 
site one hole occupies the highest e^-orbital, ~ — z^), 
i.e., the occupied orbitals (~ — z'^) or ~ \y'^ — z^)) 
alternate in all directions. This ordering and the crystal 
structure are shown in Fig. [TJ 

As a method for studying the electronic structure of 
KCuFs and the super-exchange mechanism we adopt the 
LDA-hDMFT approach [16]. Following the procedure 
presented in Ref. [3], we first calculate the LDA band- 
structure using the A/'^^-order mufhn-tin orbital method 
(NMTO). We find filled 0-bands divided by a gap of 
- 0.8 eV from the d-bands (Fig. [2]) of width Wd 3.2 eV 
- 1 eV, We^ - 2.9 eV). The energies of the d- 
crystal- field orbitals (wrt Fermi level) are —2.01 eV, 
-1.82 eV, and -1.74 eV {t2g) and -1.39 eV and -0.34 eV 
{eg). The t2g states are completely filled, do not hybridize 
with the e^-levels and thus are likely unimportant for 
orbital-ordering [l^. For the active states we construct 
a basis of localized eg NMTO Wannier functions [3] . The 
corresponding eg Hubbard model is 

im 

H~ 2 ^ ^ {Um,m' J^a,a')^ima^im' a' •) (1) 
i'm{y^m')aa' 

where riima = ^Ima^ima ^ima crcatcs an electron with 
spin cr in a Wannier orbital |m) = \x'^ —y'^) or |3z^ — 1) at 
site i; the direct and exchange [13] terms of the screened 
on-site Coulomb interaction are Umm'= U—2J{l — 5m,m') 
and J. We solve ^ using dynamical mean- field the- 
ory in the single-site approximation (DMFT) [3] and its 
cluster extension (CDMFT) |20|, using quantum Monte 
Carlo [21] as impurity solver and working^with the full 
self-energy matrix Smm' in orbital space [3]. We obtain 
the spectral matrix on the real axis by analytical con- 
tinuation We use as parameters J = 0.9 eV and 
vary U between 7 and 9 eV. These values are closeto the 
theoretical estimates based on constrained LDA [6] . 

In the paramagnetic phase, single-site DMFT calcula- 
tions yield a Mott gap of about 2.5 eV for = 7 eV, and 
4.5 eV for U = 9 eV. The system is orbitally ordered. 
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FIG. 2: (Color online) LDA 3d-bands (in eV) of KCuFs for 
the real crystal (R) and less distorted structures (see Fig. [1]). 
Blue: e^-bands. Red: t2^-bands. The Fermi level is at zero. 



and the 00 is a- type as the distortion pattern; static 
mean-field (LDA+/7, HF) calculation BSEs^ give sim- 
ilar orbital-order, however also antiferromagnetism. We 
define the orbital polarization p as the difference in occu- 
pation between the most and least occupied natural or- 
bital (diagonalizing the eg density- matrix). It turns out 
that to a good approximation p is given by the difference 
in occupation between the highest (|2)) and the the low- 
est (|1)) energy crystal- field orbital, defined in table H 
In Fig. [3] we show p as a function of temperature. We 
find that the polarization is saturated {p ^ 1) even for 
temperatures as high as 1500 K. We obtain very similar 
results in two-site CDMFT calculations. 

Using second-order perturbation theory, we calculate 
the exchange-coupling constants for the orbitally ordered 
state found with DMFT, and obtain [24] 

. 4|t-;^p(?7 + A) _ + 2J 

^ ([/ + A)2 - J2 [/ + A - 3 J [/ + A - J ' 

where fjj, are the hopping integrals from site i to site 
z', and j,/ = 1,2 are the eg crystal-field states. As 
shown in table H the calculated exchange couplings are 
in very good agreement with experimental findings. Thus 
our method gives both the correct orbital-order and the 
correct magnetic structure. 

To understand whether this orbital-order is driven by 
the exchange coupling or merely is a consequence of the 
crystal-field splitting, we consider hypothetical lattices 
with reduced deformations, measured by 7 = c/\/2a 
(tetragonal distortion) and 5 = {I — s)/{l -\- s)/2 (Jahn- 
Teller deformation). To keep the volume of the unit cell 
at the experimental value, we renormalize all lattice vec- 
tors by (7/0.95)"^/^. We calculate the Hamiltonian for a 
number of structures reducing the distortion of the real 
crystal [13] with 7 = 0.95 and S = 4.4% to the ideal 
cubic structure 7=1 and (5 = 0. The bands for some 
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Experimental structure [13] 



Imn f{\ f{\ tl 



Ley 



001 -83 -161 -161 -343 
010 98 -352 -17 59 
100 98 -17 -352 59 



57 53[25,26] 
-4 -0.4[26J 
-4 -0.4[26] 



Ideal cubic structure /o 



7 J.*' 



J—'- (th.) J+- (th.) 



001 -93 -163 -163 -282 
010 188 -327 -1 
100 188 -1 -327 



-162 
-40 
-40 




60 
60 



Crystal- field splittings A2,i and cos^ 
R RoA Ro Iaa Ioa lo.i I0.05 lo 

A2,i 1050 216 180 845 89 46 22 11 
cos 6* 0.99 0.950 0.866 0.97 0.97 0.97 0.97 0.97 0.97 

TABLE L Hopping integrals t^.'^., in the crystal- field basis 
from a site i of type 1 to a neighboring site i of type 2 
in direction /x+ny+mz. Jse are the magnetic superexchange 
couplings for the experimental structure and some representa- 
tive orbital-superexchange couplings for the ideal cubic struc- 
ture [2^. The crystal-field states are |1) = cos6'|3/^ - 1) + 
sin6'|s2 -z^) and |2) =-sin6'|3/^ - 1) + cos 0\s'^ - z^) where 
s (/) is the direction of the short (long) bond (s = x, l = y for 
a site 1). The crystal- field splitting and cos^ are given for all 
structures. All energies are in meV. 



of these structures are shown in figure [2l We use the 
notation Rs for structures with the real tetragonal dis- 
tortion 7 = 0.95 and Is for ideal (7 = 1) structures. 
The distortions affect the hopping integrals, both along 
(001) and in the x?/-plane, as shown in table [H The main 
effect is, however, the crystal- field splitting A2,i which 
decreases almost linearly with decreasing distortion, as 
expected for a Jahn- Teller system. For each structure 
we obtain the Hamiltonian i^^DA e^-bands and 

perform LDA+DMFT calculations for decreasing tem- 
peratures. At the lowest temperatures we find a- type 
00 with full orbital polarization for all structures (see 
Fig. [3]). At 800 K, the situation is qualitatively different. 
Orbital polarization remains saturated when reducing S 
from 4.4% to 1%. For smaller distortions, however, p 
rapidly goes to zero: For S = 0.2% and 7 = 1, p is 
already reduced to ~ 0.5, and becomes negligible in the 
cubic limit. Thus super-exchange alone is not sufficiently 
strong to drive orbital-ordering at T > 800 K. 

From the temperature dependence of the orbital polar- 
ization we can determine the transition temperature Tkk 
at which the Kugel-Khomskii superexchange mechanism 
would drive orbital-ordering, and thus disentangle the su- 
perexchange from the electron-phonon coupling. For this 
we study the ideal cubic structure, introducing a negli- 
gible (1 meV) crystal- field splitting as an external- field 
to break the symmetry. We find a phase transition to an 
orbitally-ordered state at Tkk ~ 350 K. The hole orbitals 



at two neighboring sites are ~ \y'^ — z'^) and ~ — ^^), 
in agreement with the original prediction of Kugel and 
Khomskii f]\. This critical temperature is sizable, but 
significantly smaller than Too ^ 800 K [lot . 

Since the screened Coulomb repulsion U is hard to 
calculate, it is crucial to identify the range of plausi- 
ble values of U and to estimate how Tkk varies in this 
range. For the experimental structure, U ^ 5 eV yields 
a tiny gap in DMFT and U ~ 6 eV a semiconducting 
gap of 1.3 eV, while KCuFa is a good insulator [l5[. It 
seems therefore unrealistic that U is smaller than 7 eV. 
It could, however, be larger. In Fig. [3] we show the re- 
sults for = 9 eV. We find a reduction to Tkk 300 K, 
reflecting the decrease in the super-exchange coupling. 
Thus we can conclude that, within single-site DMFT, 
Tkk ~ 300 — 350 K, sizable but at least a factor two 
smaller than Too 800 K. 

The ordering temperature Tkk might be even over- 
estimated by the single-site DMFT _ap proximation, as 
is common for mean-field theories [27]. To investigate 
the effects of short-range correlations, we therefore per- 
form two-site CDMFT calculations for the cubic struc- 
ture. We use a super cell containing 8 formula units, 
with axis a'= 2x, b'= 2y, c'= 2z and a two-site clus- 
ter which averages the cubic directions, i.e., imposing 
_ j^i,i±x _ Y^i,i±y _ j^i,i±z ^ nearest neighbors i 

and i' in the supercell. For U=7 eV we find that the po- 
larization starts to increase around 300-350 K, somewhat 
below the single-site transition. 

To compare these results with super-exchange theory, 
we calculate the orbital super-exchange coupling assum- 
ing that no long-range magnetic order is present [2^. 
Two representative couplings are 



J; 



SE 



(-1) 



2/7 + 6J 



U^J U-3J 



2J 



U-J U-3J 



U-J U^J 



where we adopt the pseudospin description of the orbital 
states [1], with r = 1/2 corresponding to |3z^ — 1) and 
r = —1/2 to — y^)' These couplings, shown in table 
I, are very anisotropic, with the largest about three times 
the magnetic exchange coupling along z. Again this sug- 
gests that Tkk should be larger than Tat = 38 K, but 
certainly smaller than Too ^ 800 K. 

All this indicates that in KCuFa the driving mech- 
anism for orbital-ordering is not pure superexchange. 
Further support comes from an accurate re-analysis of 
LDA+/7 results. To do this we first perform LDA+/7 cal- 
culations for the cubic and distorted structure |7| , and ob- 
tain results in agreement with previous literature [gI. [23|. 
For the undistorted KCuFa we find a fully orbitally po- 
larized solution, and the occupation matrix only slightly 



4 



Q 

Q. 









R 












'o.05^ 













500 



T/K 



1000 



1500 



FIG. 3: (Color online) Orbital polarization p as a function 
of temperature calculated with LDA+DMFT (black, blue, 
green). Black: U — 7 eV, experimental structure. Blue: 
U = 7 eV, idealized structures Rs and Is (see Fig. [1]) with de- 
creasing crystal- field. Green: U = 9 eV, Iq only. Red: cluster 
DMFT for the experimental structure and U — 7 eV. 



differs from that obtained for the experimental structure. 
This shows that the energy gain due to the distortions, 
~ 180 meV per formula unit, cannot be ascribed 
to the orbital polarization itself, but rather is an esti- 
mate of the electron-phonon coupling, enhanced by self- 
interaction correction [7] . In order to estimate the energy 
gain due to orbital polarization, we performed several 
LDA-\-U calculations for the cubic structure with differ- 
ent (fixed) occupation matrix. We find that the energy 
gain is ~ 90 meV, only half of AE. 

In conclusion, we have calculated the Kugel-Khomskii 
transition temperature for KCuFs and find a remarkably 
large Tkk ^ 350 K. Nevertheless, the super-exchange 
mechanism is not sufficiently strong to explain Too ^ 
800 K: at such a transition both super-exchange and 
electron-phonon are of comparable importance. The as- 
signment Too ^ 800 K [10] is however based on the 
temperature dependence of the orbital peak intensity sig- 
nals measured by resonant x-rays scattering [10]. A di- 
rect measurement of the evolution of the distortions with 
temperature would be highly desirable. Should orbital- 
order persist till melting, the electron-phonon coupling 
contribution would be the dominant mechanism. 

We acknowledge discussions with and helpful com- 
ments from D.I. Khomskii and P. Ghigna. The calcu- 
lations were performed on the Jiilich BlueGene/L under 
account JIFF2200. 
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